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Phenomenology from 100 large lattice 




Tanmoy Bhattacharya and Rajan Gupta a 

a T-8 Group, MS B285, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 U. S. A. 

We present a status report on simulations being done on 32 3 x 64 lattices at (3 = 6.0 using quenched Wilson 
fermions. Phenomenologically relevant results for the spectrum, decay constants, the kaon B-parameter Bk, Br, 
Bg, semi-leptonic and B — ¥ K*f form factors are given based on a statistical sample of 100 configurations. 



1. LATTICE PARAMETERS 

The 32 3 x 64 gauge lattices were generated at 
p = 6.0 using the combination of 5 over-relaxed 
(OR) and one Metropolis or Pseudo-heatbath 
sweep. We have stored lattices every 2000 OR 
sweeps. Quark propagators, using the simple 
Wilson action, have been calculated with peri- 
odic boundary conditions in all 4 directions for 
two kinds of extended sources - Wuppertal and 
Wall - at k = 0.135 (C), 0.153 (S), 0.155 (Ui), 
0.1558 (U 2 ), and 0.1563 (U 3 ). These quarks cor- 
respond to pseudoscalar mesons of mass 2816, 
977, 687, 541 and 428 MeV respectively where 
we have used 1/a = 2.314 GeV for the lattice 
scale. The three light quarks allow us to extrap- 
olate the data to the physical isospin symmetric 
light quark mass m, while the C and S k values 
are selected to be close to the physical charm and 
strange quark masses. 

We analyze three types of hadron correlators 
distinguished by the type of source/sink used 
to generate quark propagators. These are (i) 
wall source and point sink (WL), (ii) Wuppertal 
source and point sink (SL), and Wuppertal source 
and sink (SS). The effective mass m e ff(t) con- 
verges to the asymptotic value from below for WL 
correlators and from above for SL and SS correla- 
tors in all hadron channels. At the final 2cr level 
the convergence is extremely slow and there exist 
correlated fluctuations lasting 3—10 time-slices. 
As shown in Fig. [IJ the pion signal persists far 



* Based on talks presented by Tanmoy Bhattacharya and 
Rajan Gupta. These calculations have been done on the 
CM5 at LANL as part of the DOE HPCC Grand Challenge 
program, and at NCSA under a Metacenter allocation. 



Figure 1. Comparison of the convergence of 
m e ff(t) for WL and SL pion (U\Ui) correlators. 



0.31 



0.3 



0.29 



T 



M eff (p = 0) versus t 



* SL 



^ $ x v> - 



x WL 



10 



20 



30 



enough to confirm that WL and SL give the same 
mass. In all cases where we can make consistency 
checks, we find that the mean of the WL and SL 
mass is a very good estimate for the asymptotic 
value. These estimates are presented in Table |l|. 

Errors are calculated using a single elimination 
Jackknife procedure with all necessary fits and 
extrapolations done independently for each ,TK 
sample. Possible artifacts due to the quenched 
approximation are not included in the analysis 
but are discussed in the review talk by Gupta [||. 

To extrapolate the data to the physical quark 
masses, we verify that m 2 K , vn\, m p , and 
are linear in the quark mass m„ defined as 
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0.5(1/k— 1 /k c ) or determined non-perturbatively 
as discussed in Ref. 0. The results are essentially 
identical since the two estimates of m q are linearly 
related to high accuracy. An extrapolation of the 
average of WL, SL, and SS pion masses to zero 
gives /.:,.. and of the ratio rn^fm to the physical 
value 0.03182 gives Ki corresponding to m: 

k c = 0.15714(1), ki = 0.15705(1). (1) 

The lattice scale, as determined from m p , is 

a -1 (m p ) = 2.314(74) GeV. (2) 

We determine K a in three ways. We extrapolate 
m^/m 2 , m,K* /m p , and m^,/m p to m and then 
interpolate in the strange quark to match their 
physical value. Using the procedure described in 
|2J , m s evaluated at 2 GeV in the MS scheme is 

k s (Mk ) = 0.1550(1), m s (M K ) = 129(04) MeV 
K S (M K .) = 0.1547(3), m s (M K .) = 151(15) MeV 
Ks (M^ )= 0.1546(3), m s (M^ )= 157(13) MeV 

All results for matrix elements are presented 
using local operators. We use the Lepage- 
Mackenzie improvement scheme in estimating 
the perturbative renormalization factors ||], The 
"boosted" coupling is defined as 

1 (plaq) 



t 2 — ( q * = Tr/a) 



0.025. 



2. SPECTRUM 

The lattice dispersion relation for hadrons de- 
viates from the continuum form E 2 = p 2 + m 2 
due to discretization errors, and is not known 
a priori for bound states due to nonpertur- 
bartive effects. This can have important con- 
sequences for the calculation of matrix elements 
involving heavy quarks. We compare four sim- 
ple ansatze for the dispersion relation in Fig. ^| 
for CXJ\. The data favor the nearest neigh- 
bor symmetric difference relativistic dispersion 
relation sinh 2 (i?/2) — sin 2 (p/2) = sinh 2 (m/2), 
which implies that the "kinetic" mass ra-T 1 = 
(d 2 E/dp 2 |p=o) _1 = sinhm. A comparison of mi 
(given by the rate of exponential fall off) and this 
m-i is given in Table |l| for the 77 and p mesons, 
and the data show a significant difference (due to 
O(ma) effects) for the charm mesons. 



Figure 2. Test of four different dispersion rela- 
tions for CUi pion correlator 
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2.1. Baryons 

The signal in wall and Wuppertal baryon corre- 
lators becomes noisy at roughly the point where 
the m e ff begin to converge as exemplified in 
Fig. ||. We therefore use the average of SL and 
WL correlators as our best estimate. These 
numbers are given in Table 1 for the case of de- 
generate quarks. Extrapolating the data for the 
3 light quarks to m, gives 

m N a = 0.482(13); m A a = 0.590(30). (4) 

Comparing these results with GF11 data [|| we 
find that our value for M^/m p lies roughly 1.5er 
higher while that for M^/m p is consistent. 

The non-degenerate data is in excellent agree- 
ment with the Gell-Mann-Okubo formula 2(Mjv+ 
Ms) = 3Ma + Ms if either S or U\ are designated 
as the strange quark and any of the Ui as the light 
quark. Similarly, the SU (6) breaking mass differ- 
ence jtia — rriN is determined to be 289(42) MeV, 
in good agreement with the experimental value of 
293 MeV. On the other hand, the two indepen- 
dent splittings in the nucleon multiplet turn out 
to be significantly lower than their experimental 
values (using a~^(m p ) and m s (m^)): 

to Sjv - m N = 185(17) MeV cf. expt. 254 MeV 
m-E N - mjv = 300(27) MeV cf. expt. 379 MeV, 
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Table 1. Meson masses in lattice units 



Tfl\ 



mi 



m, 2 



cc 


1.217(1) 


1.541(2) 


1.229(1) 


1.564(2) 


cs 


0.854(1) 


0.962(2) 


0.880(1) 


0.999(2) 


CIJ X 


0.814(1) 


0.906(2) 


0.841(1) 


0.945(3) 


cu 2 


0.799(1) 


0.886(3) 


0.827(2) 


0.924(3) 


cu 3 


0.790(2) 


0.873(4) 


0.817(2) 


0.911(5) 


ss 


0.422(1) 


0.435(1) 


0.506(1) 


0.530(3) 


SUi 


0.364(1) 


0.372(1) 


0.465(2) 


0.482(4) 


su 2 


0.339(1) 


0.347(1) 


0.449(2) 


0.464(4) 


su 3 


0.324(1) 


0.331(2) 


0.439(3) 


0.453(5) 




0.297(1) 


0.302(1) 


0.422(3) 


0.435(5) 


UtU 2 


0.267(1) 


0.271(2) 


0.405(4) 


0.416(7) 




0.248(1) 


0.251(2) 


0.394(5) 


0.404(8) 


u 2 u 2 


0.234(1) 


0.237(2) 


0.387(5) 


0.396(9) 


u 2 u 3 


0.211(1) 


0.214(2) 


0.376(7) 


0.382(12) 


u 3 u 3 


0.185(1) 


0.187(2) 


0.363(9) 


0.365(15) 



Table 2. Baryon masses in lattice units 

SSS U1U1U1 u 2 u 2 u 2 u 3 u 3 u 3 ~ 

N 0.789(03) 0.641(04) 0.579(06) 0.540(12) 
A 0.835(04) 0.706(08) 0.660(13) 0.631(27) 



and a 20% lower still if m a (mpc) is used instead. 
Similarly, ma — m& = 322(31) MeV (cf. expt. 
440 MeV). 

3. DECAY CONSTANTS 

There are many different ways of combining SL 
and SS correlators to extract fps an d /vj some 
of which are described in Ref. [|l|. We find that all 
methods give essentially identical results, and the 
data for fps is independent of the momentum, a 
necessary condition for the reliable extraction of 
semi-leptonic form factors. We quote the mean 
value in Table |^ for two common renormaliza- 
tion schemes, naive (v^2«) and Lepage-Mackenzie 
tadpole improved (imp). The large dependence of 
fps and fv on the scheme for heavy quarks due to 
O(ma) corrections underscores the need for deter- 
mining the renormalization constant accurately. 
The results below are with the "imp" scheme. 

We find f„ = 141(6) MeV, roughly 2a larger 
than the physical value 130.7 MeV, in contrast to 
the result from GF11 collaboration [|[. Using / w 



Figure 3. WL and SL effective mass versus t for 
nucleon and A using XJ\XJ\U\ quarks. 
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to set the scale we find a _1 (/ T ) = 2.162(62) GeV. 
Similarly, we find that f K = 162(5) MeV 
(f K = 166(5) MeV) with m s set by M K (M^) as 
compared to the physical value 160 MeV. Some 
of the systematic errors cancel in the ratio fjc j f^, 
which we determine to be 1.15(1) (1.18(2)). 

The D meson with our heavy C quark has 
mi = 1789(56) and m 2 = 1972(62) MeV using 
a~ 1 (m p ), and about 7% lower using a _1 (/ 7r ). The 
physical value rap, = 1869 MeV lies within the 
range of these systematic uncertainties. There- 
fore, we take C to be the physical charm quark, 
and present results for fp> and fp a in Table 

The data in Fig. ^ show that fy 1 is qual- 
itatively different for degenerate versus non- 
degenerate quark masses. We extrapolate the 
4 lightest degenerate combinations to m to get 
fp 1 = 0.357(8), somewhat smaller than the ex- 
perimental value 0.398 and consistent with pre- 
vious determinations H. For mesons involving 
strange and charm quarks we get : 

f~l = 0.326(12), f- 1 = 0.317(7), 
fn] = 0.194(4) , f~l =0.167(6). 

Using m a (rriK) reduces these by less than la. 
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Table 3. f v and 1/fv in lattice units 



1// 



,•;< /•.' p. 

V 



0.125(2 
0.100(2 
0.091(2 
0.088(2 
0.086(2 
0.089(1 
0.083(1 
0.080(1 
0.079(1 
0.077(1 
0.074(1 
0.073(1 
0.071(1 
0.070(1 
0.068(1 



0.202(3 
0.132(2 
0.118(2 
0.112(2 
0.109(3 
0.096(1 
0.087(1 
0.084(1 
0.081(2 
0.079(1 
0.076(1 
0.073(1 
0.072(1 
0.070(1 
0.067(1 



0.104(2 
0.129(2 
0.125(3 
0.122(3 
0.120(4 
0.254(4 
0.268(5 
0.271(6 
0.272(6 
0.291(5 
0.298(6 
0.302(6 
0.308(7 
0.311(8 
0.312(9 



0.173(3 
0.176(3 
0.166(3 
0.161(4 
0.157(5 
0.282(5 
0.289(5 
0.290(6 
0.289(7 
0.307(5 
0.311(6 
0.310(7 
0.314(7 
0.315(8 
0.315(9 



Figure 4. fy 1 versus m q . 
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Table 4. fu and fr> a in MeVwith scale from m f 





m a (M K ) 


m a (M 4 ) 




mi m% 


m.\ m,2 


Id 


240(10) 229(10) 




fn/U 


1.70(4) 1.62(4) 




Sd. 


272(08) 258(08) 


279(06) 264(06) 


/d. s / /d 


1.13(2) 1.13(2) 


1.16(3) 1.15(3) 



of the LL operator (ignoring chiral logs) is 



4. Bjct Bt, B$ with Wilson fermions 

The methodology used to calculate the B- 
parameters is presented in ^ 0, so we only 
summarize the new results here. The calculation 
is done with the two operators •(/'7574V' and •075V' 
to create and annihilate the kaons. The conver- 
gence of the matrix elements is from opposite 
directions for these two cases and thus provides a 
check on the results. We find that the two sets of 
data are essentially identical except for U2U3 and 
U3U3 combinations for the LL operator. The final 
values are taken to-be the mean of the two cases 
and given in Table |5p. All data are obtained using 
1-loop improved operators as defined in Ref. (6J 
and evaluated with g 2 oost = 1.616 and fj,a = tt. 

Bk- The chiral behavior of the matrix element 



2 The numbers presented in Table 3 in LAT93 are wrong. 
The actual data is consistent with the present analysis. 



Oll \K°) = a + /3m 2 K + yp iPf + S im 4 K 

+ S 2 m 2 K pip f + S 3 (pip f ) 2 + . . .(.5) 



A fit to the data for the lightest 10 mass combina- 
tions is shown in Fig. ^. (Very similar values for 
the six parameters are obtained from fits to 4 or 
6 lightest mass combinations.) We find that even 
though the statistical quality of the data is very 
good (see Table ||), the three Si are not well deter- 
mined; only 82 is significantly different from zero. 
Terms proportional to a,/3 and S\ are pure lattice 
artifacts due to the bad chiral behavior induced 
by the mixing of the AS = 2 4-fermion operator 
with wrong chirality operators. The coefficients 
"f,S 2 ,S 3 contain artifacts in addition to the de- 
sired physical pieces and we do not have a way 
of resolving these two contributions. We simply 
assume that the 1-loop improved operator does 
a sufficiently good job of removing these residual 
artifacts. Using 1-loop matching between lattice 
and continuum gives Bk at 7.26 GeV to be 



B K = 7 + (S 3 + S 3 )m 2 K = 0.67(11) 



(6) 



a value consistent with the much more precise 
staggered fermion data [Q]. A second way of 
extracting Bk is to use Eq. |s] for pairs of points 
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Figure 5. Fit to the Bk data. 
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with different momentum transfers q: 

(E 1 B K (q 1 ) - E 2 B K (q 2 ))/(E 1 - E 2 ) = 

7 + S 2 m 2 + S 3 m(E 1 + E 2 ). 

We correct the term S 3 m(Ei + E 2 ) using the 
value of 8 3 extracted from the fit. The results 
of this analysis for the 10 light mass combina- 
tions are given in the third column of Table 0. 
Interpolating to nix, we get Bk = 0.68(12). 

The general form of Eq. § follows from Lorentz 
symmetry as m and pi-pf are the only invariants. 
The theoretical analysis of what terms are arti- 
facts in heavy-light mesons is not yet complete; 
indications are that all terms contribute. (The 
statement that a,(3,S\ are lattice artifacts relies 
on as applied to light-light systems only.) 

Since data for CU{ show almost no variation with 
momentum transfer, we therefore extrapolate the 
unsubtracted value to n\ and get Bp = 0.73(2). 

The values of B-j and Bg at the kaon mass are 

B 7 = 0.76(1), £8 = 0.81(1). (7) 

One loop running gives B s (2 GeV) « 0.67 due 
to the large anomalous dimension of Og. These 
estimates are significantly lower than those found 
in Ref. and the difference is due to a change 



Table 5. B K , B 7 and B s 









i>7 


i>8 




In — n) (n — 2) 


subtr. 






cc 


0.87(1) 0.88(2) 




0.93(1) 


0.96(1 


cs 


0.79(1) 0.81(2) 




0.92(1) 


0.95(1 




0.77(1) 0.78(3) 




0.91(1) 


0.95(1 


cu 2 


0.75(1) 0.76(3) 




0.91(2) 


0.94(2 


cu 3 


0.74(2) 0.74(4) 




0.90(2) 


0.94(2 


ss 


0.54(1) 0.61(2) 


0.79(35) 


0.87(1) 


0.91(1 




0.44(0) 0.55(2) 


0.76(27) 


0.84(1) 


0.89(1 


su 2 


0.38(1) 0.51(3) 


0.73(24) 


0.83(1) 


0.88(1 


su 3 


0.34(1) 0.47(3) 


0.70(22) 


0.82(1) 


0.87(1 




0.26(1) 0.45(3) 


0.72(19) 


0.82(1) 


0.86(1 


UiU 2 


0.14(1) 0.39(3) 


0.70(16) 


0.80(1) 


0.84(1 


UtU 3 


0.04(1) 0.34(4) 


0.68(15) 


0.78(1) 


0.83(1 


u 2 u 2 


-0.04(1) 0.32(5) 


0.69(14) 


0.78(1) 


0.82(1 


u 2 u 3 


-0.21(1) 0.26(6) 


0.68(12) 


0.76(1) 


0.81(1 


u 3 u 3 


-0.49(2) 0.20(8) 


0.67(11) 


0.74(1) 


0.79(1 



in lattice size and in the behavior versus ra q . 
Phenomenologically, a value smaller than 1.0 for 
£?8 raises the estimate for e'/e. 

5. SEMILEPTONIC FORM FACTORS 

Results for the form factors for D —¥ Kev, D —¥ 
itei>, D —¥ K*ev, D —¥ pev decays are given in Ta- 
bles ^ and 0. Similar results have also been 
presented bythe APE |], UKQCD and 
Wuppertal |ll| collaborations at this meeting. 
The details of our method are given in Ref. H Jj] 
and the main features of our analysis are: 

• The _D-meson is created at zero spatial mo- 
mentum and the momentum of the final state 
meson varies from to it /8a. This provides 
a large enough range in the invariant mass 
of the leptonic subsystem (— Q 2 ) to test the 
pole dominance hypotheses. To extract form- 
factors at Q 2 = 0, we make two kinds of fits 
using f(Q 2 ) = /(0)/(l - Q 2 /M 2 ). In the 
"pole" fits we use the lattice measured value 
for M, while in "best" fits both /(0) and M 
are free parameters. 

• Either LS or SS 2-point correlators can be used 
to cancel the matrix elements at the source and 
sink, and the exponential fall off in time as ex- 
plained in Ref. J^]. Our data show that the 
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Table 6. Semi-leptonic form factors at Q 2 = 
D -> K, K* 





m s 


m K ) 


m s (m^) 




pole 


best 


pole 


best 


u 


0.79(03) 


0.71(05) 


0.79(03) 


0.72(04) 


fa 


0.72(02) 


0.73(03) 


0.73(02) 


0.74(03) 


hlU 


0.93(01) 


1.02(04) 


0.94(01) 


1.02(03) 


V 


1.33(07) 


1.33(10) 


1.35(07) 


1.34(09) 


A a 


0.86(04) 


0.86(04) 


0.86(04) 


0.87(04) 


Ax 


0.68(03) 


0.75(05) 


0.70(03) 


0.77(04) 


A 2 


0.55(15) 


0.64(19) 


0.58(14) 


0.66(18) 


A 3 


0.86(05) 


0.86(05) 


0.87(04) 


0.87(04) 


V/A 1 


1.82(08) 


1.75(10) 


1.81(07) 


1.75(09) 


Ao/Ai 


1.15(07) 


1.15(08) 


1.14(07) 


1.13(07) 


Aij Ax 


0.76(17) 


0.85(23) 


0.78(15) 


0.87(21) 


A 3 /A x 


1.13(08) 


1.15(09) 


1.11(07) 


1.13(08) 


Table 7. 


Semi-leptonic form factors 


at Q 2 = 




D -> ir,p 


D s -> <f> 






pole 


best 






/+ 


0.70(04) 


0.56(08) 






h 


0.61(03) 


0.62(05) 






hlU 


0.91(02) 


1.07(09) 






V 


1.20(11) 


1.18(15) 


1.36(4) 




A Q 


0.83(06) 


0.82(06) 


0.80(2) 




A! 


0.59(03) 


0.67(07) 


0.67(1) 




A 2 


0.39(21) 


0.44(24) 


0.47(6) 




A 3 


0.81(06) 


0.81(06) 


0.79(2) 




V/A t 


1.89(13) 


1.77(16) 


2.01(4) 






1.29(12) 


1.23(14) 


1.19(3) 




Ail A x 


0.65(26) 


0.67(31) 


0.68(7) 




As/Ax 


1.23(14) 


1.22(14) 


1.15(3) 





meson energies obtained from LS and SS corre- 
lators agree in all cases except p = (1,1,0) and 
(2, 0, 0) pions. In these cases the 2er difference 
in E (or equivalently the amplitude) translates 
into a 2er difference in / + and /o. In the vector 
decay channels the agreement is much better. 

• The Q 2 = point lies between momentum 
transfers p = (1,1,0) and p = (2,0,0). Thus, 
the estimates of the form-factors at Q 2 = are 
sensitive to data at large momentum transfers. 

• /_)_, Ax and A 2 at Q 2 = show a significant 
difference between "pole" and "best" fits. The 
data prefer values of M which are smaller than 
that measured from 2-point correlators. Since 



?t»2 > mi, using m 2 does not reduce the differ- 
ence; it increases both estimates by about 5%. 
We take the "best" fit as our preferred value. 

• The p dependence in Zjl and Zy (O(pa) ef- 
fects) could account for the difference between 
"pole" and "best" fits. Therefore, the Lepage- 
Mackenzie [j3| tadpole improved Z's need to be 
extended to non-zero momentum transfers. 

• In order to consistently take into account all 
O(ma) effects one needs to derive the relation- 
ship between matrix elements and form-factors 
using the lattice symmetries. We have not done 
that and therefore present all results based on 
the continuum expressions. A check of this 
assumption is that the large t behavior of the 
ratio of the 2-point functions, (ViV{) / (VjVj) 
where Vi, Vj are the spatial components of the 
local vector current, is (m 2 +p 2 )/(m 2 + p 2 ). 
The data verify this at the 2% level for all mass 
combinations and momenta. 

• Our data shows that A (0) « A 3 (0) as pre- 
dicted by Lorentz invariance. 

• In most cases the statistical errors are larger 
than the difference due to the choice of m s . 
Our preferred values are with m s from m^. 

• We can use the CUx — > U\U\ measurement to 
estimate the D s —¥ <f> form-factors, assuming 
that the Zweig suppressed hairpin diagram can 
be ignored. (Ux mass is very close to m a (mpc).) 
The results given in Table show that the D — ¥ 
K* and D s —± <f> form factors are very similar. 

To summarize, our estimates of form-factors 
for charm decays are in good agreement with ex- 
perimental data. Systematic errors due to using 
the quenched approximation and due to 0(a), 
O(pa), and 0(m q a) errors remain to be under- 
stood and quantified. To extend the analysis 
to the phenomenologically interesting case of B 
mesons we propose to combine NRQCD b quarks 
with existing light propagators. 

6. THE RARE DECAY B -> K*j 

This decay is characterized by three form fac- 
tors, Ti, T 2 , and T 3 . Of these only T a (0) = T 2 (0) 
can be measured in experiments. We use 3 
methods to extrapolate our C quark data to raj 
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Figure 6. A typical example of "pole" versus 
"best" fits to /+ data. 
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using scaling relations from heavy quark effec- 
tive theory (HQET), and to Q 2 = using pole 
dominance. 

• Method I: T 2 can be measured at zero 
recoil, so we extrapolate the data to the 
physical B-meson using the HQET relation 

1 /2 

T2 (fixed recoil) nig = constant. We then 
use pole dominance with the measured value 
m* B = 5.74 GeV to estimate T 2 (0) = O.f f (1). 

• Method II: HQET along with the assump- 
tion of pole dominance predicts T 2 (0)m£ / ' 2 = 
constant. Therefore we first extrapolate T 2 to 
Q 2 = using pole dominance at the charm 
mass and then using this scaling relation to 
estimate T 2 (0) = 0.09(f) at the B mass. 

• Method III: HQET and the pole dominance 
hypothesis also predicts 7\ (0)m£ 2 = constant. 
The procedure of Method II then gives 
Ti(0) = 0.25(2) at the B mass. 

Methods I and II agree, whereas Method III 
gives a much larger value. The contradiction 
between Methods II and III arises entirely from 
the extrapolation to b mass; the lattice values 
Ti(0) = 0.39(2) and T 2 (0) = 0.4f (3) are in agree- 
ment with the theoretical prediction, based on 
Lorentz invariance, that Ti (0) = T 2 (0) for all 



Figure 7. Comparison of T\ and T 2 data with 
pole dominance expectation. 




* 1/T a - 1.0 
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l-Q 2 /mg. 

values of m q . The problem could be either in the 
pole dominance hypotheses or in using the lowest 
order HQET relations. The data in Figure show 
that Ti prefers a smaller value for to 2 -,, while T 2 
a larger one. The difference between "pole" and 
"best" fits, while small at c mass, could be large 
at the b mass. The corrections to HQET can be 
large considering the factor of 2 — 4 change in T 
due to the extrapolation. To summarize, these 
possibilities need to be brought under control be- 
fore reliable results for B —¥ K*^ can be obtained. 
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